One-loop electroweak corrections for polarized M0ller scattering 
at different renormalization schemes and conditions 
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Using two different approaches, we perform updated and detailed calculations of the complete one- 
loop (Next-to-Leading Order (NLO)) set of electroweak radiative corrections to the parity violating 
e~e~ — > e~e~ (7) scattering asymmetry. Our first approach, more classical, relies on calculations "by 
hand" with reasonable approximations. Our second approach relies on program packages FeynArts, 
FormCalc, LoopTools, and FORM. The detailed numerical analysis of the various contributions is 
provided for a wide range of energies relevant for the ultra-precise 11 GeV MOLLER experiment 
planned at the Jefferson Laboratory, as well as future experiments at the International Linear 
Collider (ILC). The numerical results obtained within the on-shell renormalization scheme using 
two different sets of renormalization conditions are in excellent agreement. We also calculate the 
total NLO correction in the Constrained Differential Renormalization (CDR) scheme. Analysis 
of the results, along with the increasing experimental precision, shows that it is feasible that the 
corrections at the Next-to-Next-to-Leading Order (NNLO) level may be important for the next 
generation of high-precision experiments. 
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I. INTRODUCTION 



The M0ller scattering measurements are not only 
one of the oldest and the best-established tools of 
modern physics, but also a clean, powerful probe of 
New Physics (NP) effects More recently, the sig- 
nificant interest to this process from both theoreti- 
cal and experimental communities has been re-ignited 
by two precision experiments: E-158 [|J at SLAC, 
which made the first observation of parity violation 
in electron-electron scattering, and MOLLER exper- 
iment planned at JLab Q. Both are dedicated to 
measuring the Parity- Violating (PV) asymmetry in 
the e~ er — > e~e~("f) scattering at low energies. The 



'Electronic address: aaleksejevs@swgc.mun.ca 



^Electronic address: svetlana. barkanova@acadiau.ca 
t Electronic address: ily@hep.by 



^Electronic address: yury@physics.berkeley.edu 
^Electronic address: vladimir.zykunov@cern.ch 



MOLLER experiment aims to measure the PV asym- 
metry in the scattering of longitudinally polarized 
electrons off unpolarized electrons with a combined 
statistical and systematic uncertainty of 36 parts per 
billion Q. With the estimated systematic contribu- 
tion of 1%, the measurement will still be statistics- 
limited, and further improvements in precision might 
be possible with additional running time. The mea- 
surement of the electron's weak charge Q^y with a 
2.3% accuracy planned by MOLLER would yield the 
most precise single measurement of the weak mixing 
angle sin 2 9w, with a fractional accuracy of 0.1%, at 
an average momentum transfer Q 2 = 0.0056 GeV 2 . 
At this precision, MOLLER can shed light on the dis- 
crepancy between the hadronic and leptonic determi- 
nations of sin 2 9w at the Z-boson pole. Furthermore, 
the difference between the values of sin 2 9w deter- 
mined at the Z-boson pole and at low Q 2 is sensitive 
to the NP effects at TeV scales. 

Before physics of interest can be extracted from the 
experimental data, radiative effects must be carefully 
treated. MOLLER's stated precision goal is signif- 
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icantly more ambitious than that of its predecessor 
E-158, so very precise theoretical input for this mea- 
surement will be crucial. In spite of significant ear- 
lier theoretical effort dedicated to calculations of Elec- 
troweak Radiative Corrections (EWC) (see the early 
review papers [H and [U, more recent Q and [7|, and 
numerous additional references in our paper |8|), we 
believe that a new level of accuracy is required for 
the next-generation, high-precision experiments. To 
match the expected experimental systematic uncer- 
tainty, it is desirable to keep the theoretical uncer- 
tainty due to the radiative corrections at or below the 
0.1% level. Obviously, calculating large sets of one- 
loop Feynman diagrams by hand is a tedious task. Re- 
cently, program packages such as FeynArts Q , Form- 
Calc LoopTools and FORM [TO] have created 
the possibility of handling the substantial number of 
diagrams reasonably quickly, minimizing probability 
of human errors, and preventing the rapid error ac- 
cumulation often unavoidable with purely numerical 
methods. One of the key features of the presented 
work is to compare, step by step, the complete one- 
loop set of EWC to the PV Moller scattering asymme- 
try calculated first by hand and then with FeynArts, 
FormCalc and LoopTools as base languages using two 
different renormalization conditions. 

FeynArts is a Mathematica package which provides 
the generation and visualization of Feynman diagrams 
and amplitudes involving Standard Model particles. 
FormCalc, a Mathematica package, reads diagrams 
generated with FeynArts and evaluates amplitudes 
with the help of the program FORM in analytical 
form. FORM, a successor to SCHOONSCHIP [U], 
is a Symbolic Manipulation System which is also es- 
sential for our computer algebra-based method, and 
is used by FormCalc as a core program. LoopTools 
provides the many-point tensor coefficient functions 
and is used to numerically evaluate scalar and ten- 
sor one-loop integrals. In FormCalc, it is possible 
for the regularization to be done either by dimen- 
sional reduction or by the usual dimensional regular- 
ization scheme. After that, one may implement one 
of the two renormalization schemes (RS), the on-shell 
scheme or the Constrained Differential Renormaliza- 
tion (CDR)scheme. For calculations done at the one- 
loop level, the CDR scheme is equivalent to regulariza- 
tion done by dimensional reduction in the MS scheme 
with redefined scale log M 2 = log /i 2 + 2, where M and 
\x are the renormalization scales in the CDR and in the 
dimensional regularization method, respectively (Io| . 

A complete automatization would limit the range 
of applications, so these packages are not "black 
box"; they require considerable human input on many 
stages. We call our approach based on FeynArts, 
FormCalc, LoopTools, and FORM "semi-automated". 
On the other hand, these packages allow modifications 
to better suit specific projects. In [l3| . for example, 



we adopted FeynArts and FormCalc for the NLO cal- 
culations of the differential cross section in electron- 
nucleon scattering. In general, the results obtained 
with these packages can be presented in both analyt- 
ical and numerical form. However, our equations for 
the EWC to the scattering asymmetry obtained with 
FeynArts and FormCalc are too lengthy and cumber- 
some, so we present only approximate equations ob- 
tained by hand. However, as we show in the numerical 
analysis section, the agreement between numerical re- 
sults obtained with the two methods - "by hand" and 
semi-automatic - is excellent. 

An additional way to ensure that our NLO EWC 
calculations are perfectly correct is the detailed com- 
parison of results calculated with different renormal- 
ization conditions within the same scheme. Of course, 
the sum of all radiative corrections forming a full 
gauge-invariant set must be independent on the choice 
of renormalization conditions. Paper fl4j , for exam- 
ple, clearly demonstrated the cancellation of gauge de- 
pendencies in one-loop corrections from self-energies, 
tadpoles, vertex and box diagrams to physical ampli- 
tudes for four-fermion processes. However, the agree- 
ment between the results evaluated in different renor- 
malization schemes can be guaranteed only if we take 
into account all orders of perturbation expansion, not 
just NLO. Since in this article we are only dealing with 
one-loop corrections, we do not expect the results pro- 
duced in different schemes to be identical. In fact, the 
difference we see between the results obtained with the 
on-shell and CDR schemes indicates the need to con- 
sider higher-order corrections. A detailed discussion 
was given in [l5[, for example. The NNLO corrections 
will be our next task. 

For now, we concentrate on achieving the best ac- 
curacy possible in one selected scheme and perform 
NLO calculations using two methods and two sets 
of renormalization conditions. For that, we choose 
the on-shell renormalization scheme with two differ- 
ent sets of renormalization conditions, the approach 
proposed by W. Hollik in (see also [13) and the 
approach suggested by A. Denner in [l8| . For brevity, 
we will call Hollik's renormalization conditions HRC 
and Denner's conditions DRC. It is obvious that any 
renormalization scheme is required to meet physical 
conditions, although it is possible to vary renormal- 
ization conditions for the sake of simplicity of the 
problem and still keep the final gauge invariant re- 
sults unchanged. As a result, contributions to the 
cross section coming from the different non-gauge in- 
variant loop corrections (self-energies and vertex cor- 
rections) could vary greatly depending on the choice 
of renormalization conditions. 

The first goal of this paper is to calculate the full 
set of one- loop EWC, both numerically with no sim- 
plifications using semi-automatic approach, and "by 
hand", analytically in a compact asymptotic form |l9l — 
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l2ll |, and compare the results. Our second goal is to 
present comparison and analysis of the various con- 
tributions to the cross-section asymmetry calculated 
within the HRC and DRC renormalization conditions. 
Our third aim is to calculate the total NLO corrections 
in the CDR scheme and estimate the importance of 
the NNLO corrections for such high-precision experi- 
ments as MOLLER. 

The rest of the paper is organized as follows. In 
Section [TT] we provide details of the basic notation, 
the lowest-order (Born or Leading Order (LO)) and 
NLO contributions to M0ller scattering. The same 
section gives a short description of photon emission 
which is essential for removal of nonphysical parame- 
ters from regularized infrared divergent cross section. 
The details of the HRC and DRC renormalization con- 
ditions and a discussion of gauge invariance can be 
found in Section Mil Analysis of analytical and nu- 
merical results in the on-shell RS using HRC and DRC 
renormalization conditions is given in the beginning of 
Section HVl Later, in the same section, the CDR re- 
sults are discussed. Section IVl includes the analysis of 
possible effects of an additional new-physics massive 
neutral boson on the observable asymmetry. Our con- 
clusions and future plans are discussed in Section IVT1 




(a) (b) 



FIG. 1: Diagrams describing nonradiative M0ller scatter- 
ing in the (a) t- and (b) u-channels. 

which depends on the Z-boson mass mz or on the 
photon mass m 7 = A. The photon mass is set to 
zero everywhere with the exception of specially indi- 
cated cases where the photon mass is taken to be an 
infinitesimal parameter that regularizes an infrared di- 
vergence. In addition, we use the functions 

\ i,k \ i,k \ i,k , \ i.k \ i.k /r\ 

A± —Mb A lT ± A 2b Mr i V 3 ! 

which are combinations of coupling constants and the 
degrees of polarizations Pb(t) of the electrons with 
momentum k\ {p\ ) given by 

M B{T) - Ay -p B (T)A A J A 2B(T) - A A -PB(T)*v J 



II. DEFINITIONS AND FRAMEWORK 

In the Standard Model, the Born cross section 
for M0ller scattering with the longitudinally-polarized 
electrons 

e-(fci) + e - (pi)->e-(A:a) + e-(p2) (1) 
can be represented in the form 

2 

= — [a 1 1 3 {u 2 D u D^ +t 2 D m D ]u ) 

+X l ] r j s 2 {D it + D iu )(D jt + D ju )], (2) 

where a = da j d cos 8 and 8 is the scattering angle 
of the detected electron with momentum k 2 in the 
center of mass system of the initial electrons. The set 
of momenta of initial (k\ and pi) and final (fc 2 and 
P2) electrons (see Fig. [T]) generates the standard set 
of Mandelstam variables, 

s = (k 1 + Pl ) 2 , t = {k 1 -k 2 ) 2 , u = {k 2 - Pl ) 2 . (3) 

We neglect the electron mass m whenever possible and 
in particular when m 2 -C s, — t, —u. 

A useful structure we employ in this paper is 

B*=—L-^ (i = 7> Z; r = *,«), (4) 
r — mr 



X v 3 = vV + a l a\ \\ j = v l aJ + oV. (6) 

Here, vector and axial-vector parts of the couplings 
have the following structure 

v t = 1, a 1 = 0, 
v z = (II + 2s 2 w )/(2s w c w ), a z - ll/{2 Sw c w ). (7) 

It should be recalled that 1% = —1/2 and sw {cw) are 
the sine (cosine) of the Weinberg mixing angle which 
is defined in terms of mz and mw according to the 
rules of the Standard Model (SM): cw = m z/ m w 
and sw = \/l — c w . The electron degrees of polar- 
ization pb(t) are labeled such that the subscripts L 
and R correspond to the values of Pb(t) — — 1 and 
Pb(t) = +1 respectively. Here, the first subscript in- 
dicates the degree of polarization for the momentum 
ki , while the second indicates the degree of polariza- 
tion for the momentum p\. Combining the degrees 
of electron beam polarizations, we can obtain four 
measurable cross sections. However by the virtue of 
the rotational invariance, two of them are identical: 
clr = ctrl- The three polarization cross sections can 
be used to construct three independent asymmetries 
(22| . Of particular interest to us is the parity- violating 
asymmetry A^r which is defined as follows 

_ OLL + OLR — CTRL — CFRR _ &LL ~ &RR 

&LL + &LR + CTRL + &RR &LL + 2&LR + O RR 

This single-polarization asymmetry corresponding to 
the scattering of longitudinally polarized electrons on 
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unpolarized electrons is proportional to the combina- 
tion 1 — 4s 2 v , and is therefore highly sensitive to small 
changes in sw That is why the asymmetry Alr was 
used as the observable in E-158 and will be measured 
in the future MOLLER experiment. At low energies 
and at Born level, the PV asymmetry A° LR is given by 



A 



2/(1 - y) 



y = -t/s.(9) 

s w 



The contribution of virtual particles (V- 
contribution) to the cross section of M0ller scattering 
is described by the three classes of diagrams: boson 
self-energies (BSE) (they include 77, jZ and ZZ 
self-energies and are shown symbolically in Fig.[2ja)), 
vertex functions (Fig. [2Kb) and Htc)), and two-boson 
exchange diagrams (boxes) shown in Fig. [3Jd, e). In 
the on-shell and CDR renormalization schemes there 
is no contribution from the electron self-energies. The 
corresponding cross section is given by the sum 



BSE 



Vcr 



(10) 



The detailed expressions for all the terms in this sum 
were given in our recent paper 

Contributions coming from the vertex correction 
graphs (with a photon in the loop), as well as the 
77 and jZ boxes suffer from the well-known infrared 
divergence. Regularization of this divergence can be 
done by giving the photon a small unphysical mass 
A. Obviously, the final result should be free of un- 
physical parameters and hence such dependence has 
to be removed. That can be done if we consider ad- 
ditional contributions associated with photon emis- 
sion diagrams (bremsstrahlung). The detailed de- 
scription of this contribution is also given in Q. The 
bremsstrahlung cross section can be broken down into 
two parts (soft and hard) as 



'IR 



•H 



(11) 



by separating the integration domain according to 
ko < u> and fco > u>, where fco is the photon energy 
(in the reference frame co-moving with the center of 
mass of the primary electrons). The parameter uj cor- 
responds to the maximum of the emitted soft-photon 
energy. First, we follow the methods ofpaper|23| to 
get a well-known result (see also |19l - l2ll . uA. l2a |) for 
the soft-photon cross section (where e is the base of 
the natural logarithm): 



°m = -(4 log — log 



A 



tu 
em 2 s 

J2 



+ 1 



1 2 J 
-log — j 

iog 2 !K. (12) 



contributions, 



a c = aY R + oft = - (4 log -= lo 



2uj , tu 



los 



7r y/s em's 

- + l-^+log 2 V. (13) 
i z 3 t 



Next, we sum the IR-terms of V- and R- 



and get a result which is free of regularization param- 
eter A. 

At this point let us continue with the discussion of 
the details of renormalization conditions we use in our 
calculations. 



III. RENORMALIZATION CONDITIONS 
AND GAUGE INVARIANCE 

To obtain the ultraviolet-finite result and render 
the parameters of the Standard Model real, we have 
to apply a renormalization procedure. For a gauge- 
invariant set, physical results should be invariant un- 
der different renormalization conditions. That is, al- 
though the contributions of the different types of dia- 
grams can vary strongly for different renormalization 
conditions, the total impact of all one-loop virtual ef- 
fects on observable quantities must remain indepen- 
dent. In other words, the contributions of separate 
self-energies and vertex correction functions strongly 
depend on the details of the renormalization con- 
ditions, and to properly account of the EWC they 
should be taken as one gauge- independent set. We 
will illustrate this for the case of the observable Aj^r, 
which is especially sensitive to the renormalization 
conditions. In addition, we can verify that our results 
are correct by comparing the computer-based (DRC) 
and "by hand" (HRC) calculations. We now briefly 
describe our two chosen renormalization conditions, 
DRC and HRC, within the on-shell renormalization 
scheme. 

Both use multiplicative renormalization constants, 
and as a result the electroweak Lagrangian, origi- 
nally written in terms of bare parameters, is sepa- 
rated into a basic Lagrangian and a counterterm La- 
grangian. The basic Lagrangian has the same form as 
the bare one, but depends on renormalized parame- 
ters and fields. The counterterm Lagrangian depends 
on renormalization constants of masses, charges and 
fields. Renormalization constants are fixed by the 
renormalization conditions, which are separated into 
two classes: the first determines the renormalization 
of the parameters, and the second fixes the renormal- 
ization of fields. The first class is related to physical 
observables at a given order of perturbation theory, 
and the second one is related to the Green's functions 
and has no effect on calculations of S-matrix elements. 
Both approaches use essentially the same renormaliza- 
tion conditions to fix the parameters of the SM in the 




FIG. 2: One-loop t-channel diagrams for the M0ller process. The circles represent the contributions of self-energies and 
vertex functions. The u-channel diagrams are obtained via the interchange ki Hp2. 



following way: 



constants: 



Retf{m 2 w ) = ReEf (m|) = Re±f(m 2 f ) = 0, 



0,p 2 



(14) 



Here, Ret^,' W (m 2 zw ) 



and T e ^ (k 2 = 0, p 2 = m 2 ) 
are the real parts of the truncated, transverse renor- 
malized boson self-energy and electron vertex correc- 
tion graphs, respectively. The longitudinal parts of 
the boson self-energy make very small contributions 
and are not considered here. The first condition of 
Eq. (|14p fixes the mass renormalization of the W-, 
Z-bosons and fermions without quark mixing. The 
second condition fixes the renormalization of electric 
charge, and is derived from the Thomson limit when 
momentum transfer k 2 = and external electrons are 
on their mass shell. As for the renormalization con- 
ditions of the fields, both approaches are quite dif- 
ferent. In HRC, field renormalization constants are 
determined from the following conditions: 



El Z (0) = 0, gf^Ej (0) = 



(15) 



However, in the DRC renormalization conditions, the 
field renormalization is defined on-shell, as it was done 
for renormalization of the SM parameters. This ex- 
plicitly introduces an additional set of conditions, be- 
sides Eq. (fT4")) and Eq. (fT5"|). which read: 



ReE 



7 z 



dk 2 



Re 



(16) 



As a result, in DRC, renormalization constants for 
the fields of vector bosons are calculated in a rela- 
tively simple way, without the mass-renormalization 



SZ, 
SZ 

sz 



(D) 

w 

(D) 

z 

(D) 
Z-y 



Re^E- (m 2 w ) 



Re 



dk 2 



l z) i 



ReEf (0) , SZ\ U Z 



(D) 



2 R.g5j/ji 



6ZW 



(17) 



They can be presented through truncated and non- 
renormalized self-energy graphs. In comparison 
with HRC, where the renormalization conditions of 
Eq. Q16[) are not present, field renormalization con- 
stants are defined in a different way and depend on 
the mass-renormalization constants: 



5ZW 



SZ 



+2 



(H) 



-W 



nr 



E^(0) 



swcw 



SZ 



H) 
W 



s w 



dk 2 

( Sm 2 z 



Sm 2 z 



:EJ(0) 



Sm 2 w 



Sm 2 w 



^z 



(0) 



Sw 



Mr 



cwsw 



SZ 



(H) 



6ZW 



(18) 



The presence of the mass renormalization constants 
in the field-renormalization Eq. (|18[) increases the val- 
ues of the truncated and renormalized self-energy di- 
agrams, and the dominant NLO contributions to the 
observable cross section come from these diagrams. 
In DRC, the mass renormalization constants appear 
in renormalization constants of the electroweak cou- 
plings, and hence we observe comparable contribu- 
tions coming from both self-energies and vertex cor- 
rections. Of course, such a comparison has no physical 
meaning since neither self-energies nor vertex correc- 
tions represent a gauge-invariant set on their own. As 
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is well known, only the sum of both groups is gauge 
invariant; later we show that both approaches give ex- 
actly the same results for the observable asymmetry. 
We would like to highlight that it is important to ex- 
ercise caution when comparing separate contributions 
arising from the different renormalization conditions. 
This point is illustrated by Fig. [31 [Hand El where one 
can see various renormalized vector boson self-energies 
calculated with both DRC and HRC. 

Let us note that the same renormalization condi- 
tions are imposed on the electromagnetic field. As 
we can see from Fig. [HJ which shows the truncated 
and renormalized 77 self-energies, there is no differ- 
ence whatsoever between the two sets of conditions. 
The situation is quite different if we look at the re- 
sults for the truncated and renormalized ZZ , 7Z and 
WW self-energies. In Fig. @] and [SJ we can see a 
substantial difference in the results obtained within 
the two sets of renormalization conditions, where the 
DRC set systematically leads to the self-energies being 
smaller in magnitude. As a result, with DRC, the self- 
energy contributions to M0ller asymmetry are roughly 
a factor of two smaller in value compared to the val- 
ues given by HRC. However, adding the DRC vertex 
corrections restores the total correction obtained with 
DRC to within 0.001% of the HRC result at all ener- 
gies relevant to the planned MOLLER experiment at 
JLab. 

Now we are ready to proceed to the analysis of the 
results. Section TlVAl is done in a manner independent 
of the renormalization conditions, and all notations 
that follow below can be applied to both DRC and 
HRC of the on-shell scheme. 



IV. RESULTS AND ANALYSIS 



case. We checked that variations of the light quark 
masses around the outlined values have a negligible 
impact on the values of the polarization asymmetry, 
so the choice of quark masses does not introduce a sig- 
nificant uncertainty to our results. An earlier work by 
@, relevant to observables measured at very low mo- 
mentum transfers, which determined the weak mixing 
angle in the MS-scheme, also argued that the uncer- 
tainty from non-perturbative hadronic contributions 
is small compared to anticipated experimental uncer- 
tainties. Although [25[ argued that the most signif- 
icant source of theoretical uncertainty on scattering 
asymmetry comes from the hadronic contributions to 
the "fZ vacuum polarization, we find that in our cal- 
culations hadronic contributions are under good con- 
trol. Finally, for the mass of the Higgs boson, we take 
uiH = 115 GeV. Although this mass is still to be de- 
termined experimentally, the dependence of the EWC 
on run is rather weak. 

Let us determine the physical impact of this contri- 
bution to the observable Alr, by defining the relative 
corrections to the Born asymmetry as 

$a = {A LR — a lr)/ a lri 

where the index C stands for a specific contribution, 
for example C = BSE, Ver, Box. Let indicies 7SE, 
7ZSE and ZSE denote 77-, 7Z-, and ZZ-BSE contri- 
butions, respectively. The main subject of the follow- 
ing analysis is "weak 9 relative corrections, which are 
defined as all BSE contributions (including the 77-SE 
which is not weak by nature, but is needed here to ac- 
count for all IR-finite contributions to the asymmetry) 
plus heavy vertex (HV) contributions ("heavy" means 
"massive", i.e. Z- or W-boson), ZZ- and H^M^-boxes. 
In summary: weak = BSE+EV+ ZZ+WW. 



A. On-Shell Renormalization 



Analysis of BSE contributions to PV asymmetry 



For the numerical analysis, we use a = 
1/137.035999, m w = 80.399 GeV, and m z = 
91.1876 GeV according to [26[, take electron, muon, 
and r-lepton masses as m e — 0.510998910 MeV, m M = 
0.105658367 GeV, m T = 1.77684 GeV, and quark 
masses for loop contributions as m u = 0.06983 GeV, 
m c = 1.2 GeV, m t = 174 GeV, m d = 0.06984 GeV, 
m s = 0.15 GeV, and m& = 4.6 GeV. The light quark 
masses provide Aa^ d (m 2 z )=0. 02757 [13, where 



A«i 5 i( S ) 



a 



Qf[ log 



(19) 



Qf is the electric charge of fermion / in proton charge 
units q, (q = V Ana). We can see that the use of 
the light quark masses as parameters regulated by the 
hadronic vacuum polarization is a better choice in this 



We start with the ZSE-contribution, where 

AZSE aQ )|ll-rr _ g |ll-rr 
rZSE _ A LR ~ A LR _ "aa+<'aa E CT o'o 



A LR 

v |ll-rr 
c°|ll-rr 



ct°|ll-rr 




(20) 



The operation £/|ll-rr under expression E means 
Ell — E RR . Subscript 00 denotes the unpolarized 
cross section. The approximate equality is possible 
because <7(x) SE /o'oo is very small. The denominator of 
the last fraction is calculated directly from Eq. @: 

<t°|ll-rr = l^a 2 v z a z s{D^ + D ru )(D zt + D Zu ) 
K-i2na 2 v z a z A T {D^ t + D- 1u ). (21) 
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FIG. 3: Truncated and renormalized 77 self-energies in both sets of renormalization conditions. The right graph shows 
the low-energy domain. 



With simplifications, the numerator of Eq. (f2U|) is 
<x zse |ll-rr « ^D zzt (D^M z ^ 



i od 



-16 7 ra 2 w z a z ^-( J D'' t + ZV U ) 
ml 



Finally. 



x(E|(t) + Ef(u)). 



rZSE ^ Sy(t) + £|(u) 

~ 2m| 



(22) 



(23) 



At small r (Eq. Q ) corresponding to £j a b = 11 GeV 
and 9 = 90°, the corrections are (Tf SE (HRC) w 0.0309 
vs. 5% SE (DRC) « -0.0105. 

Similarly, for 7ZSE-contribution 

■2 



V LL-RR 



—D 1 s zt [D' lt {M2^ Zl + M^ 777 ) 



-D^{M22 Zl + A/ Z d 777 )]| LL _RR + (t o u) 
l 67ra 2 a z (^ 7 * + £ 7 ") ( + ) ( 24 ) 



and 



For 7ZSE at = 11 GeV and 9 = 90°, the correc- 
tions are <5 7ZSE (HRC) « -0.6028 vs. <^ SB (DRC) 
w —0.2909. It is important to note that the deviation 
in TiTj? has a dramatic impact on such a sensitive ob- 
servable as Alr. For example, the uncertainty in 
of 1% will result in a change in (J 7 ^ E of up to 0.05. 

All terms with properties of Eq. (|2"0"]) contribute ad- 
ditively to the total correction, for example, 



nZSE+ZSE 



c-yZSE 



(26) 



We call such contributions additive. 7SE gives a non- 
additive and small contribution that we consider later. 



2. Analysis of HV and box contribution to PV 
asymmetry 

Starting with the A2-contribution, which comes 
from the triangle diagrams with an additional mas- 
sive boson, Z or IV, we get 



jA 2 _ o- - 1 LL-RR 
c°|ll-rr 



(27) 



The numerator, with some approximations, is 

<t A2 |ll-rr « 8a 3 v z a z s(D~< t + D 1U ) 
'A 2 (t,mz) A 2 (!i,mz)- 



t 



(28) 
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FIG. 4: Truncated and renormalized ZZ and WW self-energies in both sets of renormalization conditions. 



so the correction is proportional to A2 in the following 
way: 



we find 



. am 2 z fA 2 (t, m z ) , A 2 (u,m z ) 

A ~ 



-in 



t 



(29) 



In HRC, we can simplify the result by using series 
expansion of A2 at small t: 



which gives 

A 2 a tu 23 
07T m z 6 

The numerical value obtained from above at E\ 



(31) 



lab — 

11 GeV and 6 = 90° gives 8^ w -0.0125, which is in 
agreement with the exact (semi-automatic) numerical 
calculations. 

The A3-contribution, which represents the triangle 
diagrams with a 3-boson vertex, WWj or WWZ, is 
calculated in a similar way, so we present only the 
final result: 



3am| / A 3 (t,mw) A 3 (u,mw) 



32ns 2 w 



v z a z \ 



t 



(32) 



After simplifications and series expansion of A3 at 
small t, 



tt9(1-4 S ^)- 



(34) 



Using Eq. (M} for E uh = 11 GeV and 6 = 90°, we 
obtain <5^ 3 w 0.0118. Again, this approximate value 
calculated "by hand" is in a very good agreement with 
the exact result obtaned with our second, computer- 
based approach. 

The box part is UV-finite and does not require the 
renormalization procedure. We divide the box contri- 
bution into QED (77- and 7.Z-boxes) and a heavy-box 
part (HB = ZZ + WW): 



QED 



'HB 



(35) 



The types of boxes are shown in Fig. [2jd, e). The 
IR-divergent QED-part of boxes (the first term in 
Eq. (|35| ) is described in detail both analytically and 
numerically in Q. For the purely- weak part of the 
boxes (the second term), the equations are derived in 
the low-energy approximation. 

The total weak correction to A^r includes the HB 
cross section: 



A 3 (t,m w ) 



+ 0{^r), (33) 



27m? w 



'HB 



E 



(B 



k 

ZZ 



(i-B-u), (36) 
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FIG. 5: Truncated and renormalized 7Z self-energies in both sets of renormalization conditions. The right graph shows 
the low-energy domain. 



where the expressions for B k ^ take a form 

B h zz = D kt \ Bk 5 zz + (D kt + D ku )\ Bk S 2 zz , 

Bww = D kt X Ck 5] vw + (D kt + D ku )\% k 5 2 w ^7) 

The combinations of the coupling constants are given 
in Eq. . Let us recall the coupling constants for the 
heavy boxes: 

v B = {v z f + (a z )\ a B = 2v z a z , 

v c =a c = l/(4s 2 w ). (38) 

At s, \t\, \u\ -C rn z , the corrections have a form: 





3u 2 
2m|' 


s zz = 


3s 2 
2m|' 


°ww 


2u 2 


' °ww 


s 2 


m w 


2m 2 w 



(39) 



At last, after simplification at small t, for the rela- 
tive corrections to PV asymmetry coming from heavy 
boxes we find: 

tZZ _ 3a B cWW _ « ( A r,\ 

°A 7T- v ) °A ~ . TTl A 2 \ - V 4U / 

2?r 4nsfo(l- 4s^J 

The numerical values obtained from the equations 
above at E lah = 11 GeV and 6 = 90° give 8f z « 
-0.0013 and (5]f ^ 0.0238, which is once again in 



good agreement with the exact results evaluated with 
help of the FeynArts, FormCalc, LoopTools and Form 
program packages. 



3. Numerical analysis on EWC to PV asymmetry 

In the table below, we present the contributions to 
relative weak corrections calculated using two different 
approaches. In the first approach, we use approximate 
and compact expressions derived "by hand" with the 
application of HRC. In the second, we use computer- 
based analytical (FeynArts, FormCalc, and FORM) 
and when numerical (LoopTools) calculations, with 
DRC. 

Table 1 demonstrates that the 77-SE contribution 
is small, non-additive and, as expected, is the same 
whether obtained in HRC or DRC. The 7Z-SE, ZZ- 
SE and HV contributions are rather sizeable, are all 
additive, and are different for HRC and DRC. The 
ZZ-box contribution is small, and the M^W^-box is 
dominant for the weak box correction. Both the ZZ- 
box and WW^-box are additive and their sum is in 
excellent agreement regardless the method of calcu- 
lations. The total relative weak correction is signifi- 
cant and in excellent agreement between the different 
methods. That confirms that we are dealing with a 
gauge-invariant set of graphs. The discrepancy be- 
tween the two approaches is ~ 0.0001 at 9 = 90°, but 
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TABLE I: The Born asymmetry A LR and the structure of relative weak corrections to it for _Bi a b = 11 GeV at different 
0. 



9° 


20 


30 


40 


50 


60 


70 


80 


90 


A° LR , ppb 


6.63 


15.19 


27.45 


43.05 


60.69 


77.68 


90.28 


94.97 


77-SE, DRC 
77-SE, HRC 


-0.0043 
-0.0043 


-0.0049 
-0.0049 


-0.0054 
-0.0054 


-0.0058 
-0.0058 


-0.0062 
-0.0062 


-0.0064 
-0.0064 


-0.0066 
-0.0066 


-0.0067 
-0.0067 


■yZ-SE, DRC 
■yZ-SE, HRC 


-0.2919 
-0.6051 


-0.2916 
-0.6043 


-0.2914 
-0.6042 


-0.2912 
-0.6038 


-0.2911 
-0.6034 


-0.2910 
-0.6031 


-0.2909 
-0.6028 


-0.2909 
-0.6028 


ZZ-SE, DRC 
ZZ-SE, HRC 


-0.0105 
0.0309 


-0.0105 
0.0309 


-0.0105 
0.0309 


-0.0105 
0.0309 


-0.0105 
0.0309 


-0.0105 
0.0309 


-0.0105 
0.0309 


-0.0105 
0.0309 


HV, DRC 
HV, HRC 


-0.2946 
-0.0015 


-0.2633 
-0.0012 


-0.2727 
-0.0010 


-0.2703 
-0.0009 


-0.2714 
-0.0008 


-0.2712 
-0.0007 


-0.2711 
-0.0007 


-0.2710 
-0.0007 


ZZ-box, exact 
ZZ-box, approx. 


-0.0013 
-0.0013 


-0.0013 
-0.0013 


-0.0013 
-0.0013 


-0.0013 
-0.0013 


-0.0013 
-0.0013 


-0.0013 
-0.0013 


-0.0013 
-0.0013 


-0.0013 
-0.0013 


VFVF-box, exact 
iyiy-box, approx. 


0.0239 
0.0238 


0.0238 
0.0238 


0.0238 
0.0238 


0.0239 
0.0238 


0.0239 
0.0238 


0.0238 
0.0238 


0.0238 
0.0238 


0.0238 
0.0238 


total weak, DRC, exact 
total weak, HRC, approx. 


-0.5643 
-0.5526 


-0.5430 
-0.5514 


-0.5508 
-0.5511 


-0.5489 
-0.5505 


-0.5500 
-0.5500 


-0.5495 
-0.5496 


-0.5493 
-0.5493 


-0.5493 
-0.5493 




-1 2 3 

10 1 10 10 10 



Vj, GeV 

FIG. 6: The relative weak (solid line in DRC (semi- 
automated) and dotted line in HRC ("by hand")) and 
QED (dashed line) corrections to the Born asymmetry 
A° LR versus y/s at 8 = 90°. The filled circle corresponds 
to our predictions for the MOLLER experiment. 



becomes larger with decreasing 9. 

In Fig. [S] we can see the relative weak corrections 
shown by solid line for DRC (exact) and dotted line for 
HRC (approximate) . The dashed line shows the QED 
correction obtained by including soft bremsstrahlung 
to the Born asymmetry A Q LR . We can see that for low 



energy region 1 < y/s < 30 GeV the results calculated 
by the two methods are in excellent agreement. It is 
worth mentioning here that the semi-automated nu- 
merical calculations of boxes in the region of y/s <C 1 
GeV suffer from the numerical instability due to Lan- 
dau singularities. As for our approximated calcula- 
tions, we have used the small-energy approximation 
with the expansion parameters taken as r/mr zw for 
energies y/s < 30 GeV. In any case, for the 11 GeV 
relevant for the planned JLab experiment, the consis- 
tency of our calculations in both approaches is obvi- 
ous, with a difference of ~ 0.01% or less. The dotted 
line for y/s > 500 GeV on the Fig. [6] is obtained us- 
ing HRC with the help of equations from (28|, which 
used the high-energy approximation. We can see good 
a agreement between our results for the high-energy 
region y/s > 500 GeV which becomes better with en- 
ergy increase. For y/s > 50 GeV we have excellent 
agreement with the result of [24| if we use their SM 
parameters (see Q). Furthermore, the relative QED 
correction (see Fig. 8 in [24j and dashed line in Fig. E 
here) is also in good qualitative and numerical agree- 
ment. In this case, we apply the same cut on the soft 
photon emission energy as in [24| (uj/y/s = 0.05). At 
the low-energy point corresponding to the E-158 ex- 
periment, and using our set of input parameters (a, 
mw and mz) we find that 8^ eak sa —54%. If we trans- 
late our input_parameters to the set a, Gf and mz 
accordin g to |24| , we obtain good agreement with the 
result of [Hi. 



11 



B. Constrained Differential Renormalization 

The CDR (Constrained Differential Renormaliza- 
tion) scheme, which provides renormalized expressions 
for Feynman graphs preserving the Ward identities, 
was introduced at the one- loop level in |3(j. [HI ex- 
pands on [30| to introduce the techniques for one-loop 
calculations in any renormalizable theory in four di- 
mensions. The procedure has been implemented in 
FormCalc and LoopTools, which allows us to evalu- 
ate NLO EWC in CDR. Since our "scheme of choice" 
at the moment is on-shell, which is more suitable for 
calculating EWC beyond one-loop, we do not provide 
the same detailed analysis and step-by-step compari- 
son between the two methods for CDR as we do for 
on-shell. The reason we evaluate NLO EWC in CDR 
is to obtain some indication of the size of the higher- 
order effects (NNLO and beyond) to see if there is 
enough motivation to do these very involved calcula- 
tions in the future. 

In Fig. we can see the relative total correction 

gtat = ^tot _ a 0y a 

to the unpolarized cross section versus y/s at 9 = 90° 
for different RS: on-shell and CDR. In the region of 
small energies, the difference between the two schemes 
is almost constant and rather small (~ 0.01), but 
grows at yfs > mz- It is well known that in the region 
of small energies, the correction to the cross section 
is dominated by the QED contribution. However, in 
the high-energy region the weak correction becomes 
comparable to QED. Since the difference between the 
on-shell and CDR results grows substantially as the 
weak correction becomes larger, it is clear that for an 
observable such as the PV asymmetry the difference 
between the on-shell and CDR schemes will be sizeable 
for the entire spectrum of energies y/s < 2000 GeV. 
Because of that, we expect that the NNLO correction 
to the PV asymmetry may become important to PV 
precision physics in the future. 

Fig. [5] shows the relative weak (lower lines), and 
QED (upper lines) corrections to the Born asymmetry 
A° LR versus -y/s at 9 = 90°. The difference is signifi- 
cant and is growing with increasing ^fs. According to 
our calculations for Ei a b = 11 GeV, w = 0.05 y/s and 9 
= 90°, the total radiative correction to PV asymmetry 
is -69.8% with on-shell and -58.5% with CDR. The 
difference is not at all surprising. For E-158, for ex- 
ample, the one-loop weak corrections were found to be 
about -40% in the MS scheme [H and about -50% 
in the on-shell scheme (2ll |25| . 

The physical, NLO-corrected asymmetries, com- 
puted in both on-shell and CDR schemes, are com- 
pared in Fig. [5] Here, for consistency with the MS 
definition of the couplings to 0(a) [32j, we use s\ = 
sin 2 6 W (M Z ) = 0.2313 [H] in the expression of the 
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FIG. 7: The relative total corrections to the unpolarized 
cross section versus y^s at 6 = 90°. The filled circle cor- 
responds to our predictions to the MOLLER experiment. 
Solid line corresponds to CDR and dotted line to on-shell 
RS. 

Born asymmetry. We find that the predictions for the 
physical PV asymmetry, computed to the same order 
in perturbation theory in two different schemes, dif- 
fer by about 3%. The difference is an indication of 
the order of magnitude the higher-order, NNLO and 
beyond, terms. 

The [25| estimated that the higher-order corrections 
are suppressed by ~ 0.1% relative to the one-loop re- 
sult, possibly 5% in some cases, and thus are not sig- 
nificant source of uncertainty. However, we conclude 
that although the corrections at the NNLO level were 
not mandated by the previously achievable experimen- 
tal precision, they may become important for the next 
generation of experiments. 



V. EFFECT OF ADDITIONAL MASSIVE 
NEUTRAL BOSON 

Let us now add a very simple NP assumption to 
our SM calculations and show how this NP contri- 
bution affects the observable asymmetry. The reason 
we want to do it in here is to investigate if the two 
complimentary methods we used in the previous sec- 
tions, "by-hand" and semi-automated, can be applied 
in the NP domain. As we mention in the Introduc- 
tion, FeynArts, FormCalc, LoopTools, and FORM are 
not "black box" programs and can be modified for 
specific projects, including adding the NP sector. As 
was already concluded in [33[ and [Hj], the proposed 
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FIG. 8: The relative weak (lower lines) and QED (upper 
lines) corrections to the Born asymmetry A° LR versus yfs at 
9 = 90° . The filled circle corresponds to our predictions to 
the MOLLER experiment. Solid lines correspond to CDR 
and dotted lines to on-shell RS. 



MOLLER measurement could be influenced by radia- 
tive loop effects of new-physics particles. This type 
of calculation is out of scope of this paper, but we 
plan to provide a full estimation in our future publi- 
cation. For now, we assume that there is just one ad- 
ditional neutral boson (ANB), or Z'-boson, with the 
usual V — A structure of interaction with fermions, 
vector (axial) coupling constants v z (a z ) and mass 
mz> ■ From the analysis done in the previous section, 
we can clearly see that in the low-energy region where 
s, \ t\, \u\ <C m \ < m z ,, contributions are mainly sup- 
pressed by propagator factors like D z r . In this sec- 
tion, our goal is to analyze the contribution of Z'-Born 
and ZZ'-box diagrams to the observable scattering 
asymmetry for MOLLER experiment. The only sig- 
nificant contribution to the Born asymmetry comes 
from the interference terms from the Z' and photon 
diagrams. The relative correction to the Born asym- 
metry coming from Z'-boson is additive, and is given 
by 



, _ v z 'a z ' 



(41) 



According to d|, the goal of MOLLER is to measure 
the PV asymmetry to a precision of 2% (0.73 ppb). 
With this uncertainty, and assuming the identical cou- 
pling constants for Z and Z', it should be possible to 
detect ANB with a mass up to mz' = \f m 2 z / {) .02 w 
7mz- The sensitivity of MOLLER to Z' increases if 



10 c 




■fs, GeV 



10 



10 



10 



FIG. 9: The NLO-corrected asymmetries vs \fs at 9 = 
90° , computed in on-shell RS (solid line) and CDR (dotted 
line). The CDR Born asymmetry uses the MS definition 
of s| = sin 2 6 w (Mz) = 0.2313 



its parity- violating couplings are larger than those of 
Z° ', making the measurements of PV complementary 
to the direct searches at high energies. 

The one-loop diagrams including ANB give signif- 
icantly smaller contributions. As an example, let us 
consider ZZ'-box. As before, we perform our calcu- 
lations by both approximate ("by-hand") and exact 
(with FeynArts and FormCalc) methods, and get an 
excellent agreement. The expressions derived as a re- 
sult of our approximate approach are presented below. 

For the ZZ'-box contribution, the cross section can 
be expressed by the following short equation: 



r ZZ'-box 



3cr 
s 



L E 

k=7,Z 



2\f k s 2 (D 



X B ' k (D kt u 2 + D ku t 2 ) 

+ (Z O Z% (42) 



where 



L = 



1 



D 



log 



ii^Z — " L Z f TYlZ 

and the functions k are expressed through 



(43) 



B' Z Z' , Z Z' 

v = v v + a a 



v z a z '+v z 'a z . (44) 



To obtain L, we calculate the master scalar integral: 

dy 



i i i 

2, 



z dz xdx 



t(x-l) + (t-m%,)(l-z) 



(45) 
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According to the methods described in Sec- 
tion HVAl we can now calculate the relative correction 
to the observable asymmetry from the ANB contribu- 
tion (i.e. from ZZ'-box) <5^ NB as: 

= ^ 4^ L. (46) 

This correction is additive and becomes less impor- 
tant with increasing mz 1 ■ However, this suppression 
is not very dramatic due to the growing log in the 
numerator. If we take the MOLLER kinematics and 
assume that v z = v z , a z = a z , then for r m = 
mz> I mz = 1 the correction is twice the contribution 
from ZZ-box: <^ NB -0.0025465. As r m grows, 
the correction decreases: at r rn = 2 the correction is 
s anb _ -0.0011768, and for r m = 10 the correction 
is S^ NB m -0.0001185. At r m = 20, the correction 
becomes completely negligible: 6% NB « -0.0000382. 
However, the possible contributions of new-physics 
particles to the M0ller scattering deserves further at- 
tention, and we intend to continue our work in this 
direction. 

One of the simplest supersymmetric SM exten- 
sions is the Minimal Supersymmetric Standard Model 
(MSSM), and it gives a useful framework for dis- 
cussing SUSY phenomenology. For e~e~ scattering, 
MSSM contributions will arise at the one-loop or- 
der, and the large suppression of the SM weak charge 
makes the weak charge sensitive to the effects of new 
physics. According to [33^ . the loop corrections in the 
MSSM can be as large as ~ 4% for the weak charge of 
the proton and ~ 8% for the weak charge of the elec- 
tron, which is close to the current level of experimental 
and theoretical precision available for the low-energy 
studies. Obviously, before we can interpret these high- 
precision scattering experiments in terms of possible 
new physics, it is crucial to have the SM EWC under 
a very firm control. 

VI. CONCLUSION 

In the presented work, we perform detailed calcu- 
lations of the complete one-loop set of electroweak 
radiative corrections to the parity violating e~e~ — > 
e~e - (7) scattering asymmetry both at low and high 
energies using the on-shell renormalization conditions 
proposed in \Wh (see also [l?}) and the conditions sug- 
gested in fl8| |. Although contributions from the self- 
energies and vertex diagrams calculated with the two 
sets of renormalization conditions differ significantly, 
our full gauge-invariant set still guarantees that the 
total relative weak corrections are in excellent agree- 
ment for the two methods of calculation. 

Obviously, it is important to exercise caution when 
comparing separate contributions arising from the dif- 



ferent renormalization conditions unless these contri- 
butions form a gauge-invariant set (like boxes). Al- 
though this is a well-known fact in principle, it is still 
useful to demonstrate this in detail numerically for a 
specific example. We hope that our results illustrating 
the structure of relative weak corrections evaluated at 
different renormalization conditions will be of educa- 
tional value to researchers staring work in this area. 

In addition, we compare the asymptotic results ob- 
tained analytically, "by hand" (with HRC), with some 
approximations, and semi-automatically (with DRC), 
with no approximations required. As a result, we have 
a good agreement for the whole < y/s < 50 GeV 
energy region. More specifically, for the kinematics 
relevant to the 11 GeV MOLLER experiment planned 
at JLab, our agreement within two approaches for the 
complete one-loop set of electroweak radiative correc- 
tions is better than 0.1%. We found no significant 
theoretical uncertainty coming from the largest possi- 
ble source, the hadronic contributions to the vacuum 
polarization. The dependence on other uncertain in- 
put parameters, like the mass of the Higgs boson, is 
extremely weak and well below 0.1%. We conclude 
that the excellent agreement we obtained between the 
results calculated "by hand" and semi-automatically 
serves as a good illustration of opportunities offered 
by FeynArts, FormCalc, LoopTools, and FORM. 

Considering the large size of the obtained radia- 
tive effects, it is obvious that the careful procedure 
for taking into account radiative correction is essen- 
tial. Our plans include the construction of a Monte 
Carlo generator for the simulation of radiative events 
within M0ller scattering to make our work directly 
useful to the experiment. Since we are now assured of 
the reliability of our calculations, we plan to base this 
Monte Carlo on the maximum-precision results from 
our semi-automatic approach. 

Although making sure that the results obtained by 
two different approaches using two renormalization 
conditions are identical assures us that our NLO EWC 
calculations are error-free, it does not address the 
question of the size of NNLO corrections. The two- 
loop corrections are beyond the scope of this work, but 
we plan to address them in the future. One way to 
find some indication of the size of higher-order contri- 
butions is to compare physical observables computed 
to the same order in perturbation theory in differ- 
ent renormalization schemes. Our calculations in the 
on-shell and CDR schemes show that while the NLO 
terms differ by about 11%, the PV asymmetries differ 
by about 3%. At the level of precision of the future ex- 
periments such as MOLLER, higher-order corrections 
become important. 

To see if the two complimentary approaches we suc- 
cessfully used for the SM calculations can be applied 
in the NP domain, we expanded FeynArts, FormCalc, 
LoopTools, and FORM to include an additional neu- 
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tral boson (Z') : calculated the relevant correction, and 
then obtained the same result by hand. Possible other 
contributions of new-physics particles to the M0ller 
asymmetry still need to be investigated, and many of 
them can be included into the program packages men- 
tioned above. 

We believe that the future experiments at JLab and 
the ILC will mandate evaluation of the EWC beyond 
one loop. Once all the SM corrections are under con- 
trol, it is worth considering NLO corrections including 
new-physics particles, starting with the Minimal Su- 
per Symmetric Model (MSSM). The most straightfor- 
ward way to address these corrections is by employing 
the CDR scheme [3(| , [M[ because the CDR approach 
can be easily expanded to MSSM. However, whether 
the CDR scheme will be applicable in evaluating the 
EWC at the NNLO level is still an open question. 
Our preliminary plan is to address the NNLO EWC 
with the on-shell scheme first, and if the effect is sig- 
nificant, stay with the same scheme for calculating 
contributions coming from the new-physics particles. 



The simple example of the ZZ'-box we consider in 
Section [V] is gauge-invariant and is thus not affected 
by the choice of renormalization, but we have to be 
careful when choosing the scheme for our future work. 
Any suggestions from the community regarding the 
best approach to this task would be greatly appreci- 
ated. 
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